Inhomogeneous Superfluid Phases in the Unitary ^Li-^°K Mixture. 
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We have determined the complete phase diagram of the ''Li-^^K mixture within mean-field the- 
ory taking the inhomogeneities of the fermion states into account exactly. We find two different 
inhomogeneous superfiuid phases in mixtures with a majority of atoms, namely the Larkin- 
Ovchinnikov (LO^) phase with one inhomogeneous direction and a cubic phase (LO^) where three 
spatial translational symmetries are broken. We determine the transition between these two phases 
by solving the Bogoliubov-de Gennes equations in the LO^ phase. Subsequently, we calculate the 
atomic density modulation of the atoms in the LO^ phase and show that for the minority atoms it 
is sufficiently large to be visible in experiment. 
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Since the first realization of superfluidity in a two- 
component gas of fermionic atoms a large number of ex- 
citing experiments have been performed. For instance, 
with the use of a Feshbach resonance the inter-particle in- 
teractions can be tuned and hereby the crossover between 
the Bose-Einstein condensation of molecules and Cooper 
pairs was studied [ll, . By inducing spin flips in a Fermi 
gas the number of particles in different spin states can 
be changed, opening up the possibility of studying the 
influence of a population imbalance on the phase tran- 
sition to a superfluid state. This possibility has indeed 
been materialized in a two-component gas of ^Li atoms 
and the phase diagram turns out to be governed by a 
tricritical point below which the gas phase separates into 
a superfluid and normal region [3-6]. 

One of the exotic states of matter that may also be real- 
ized in a system of ultracold Fermi atoms is a superfluid 
where the Cooper pairs have a nonzero center-of-mass 
momentum. The first to propose this possibility in the 
context of superfluid films in a magnetic field were Fulde 
and Ferrel and independently Larkin and Ovchinnikov 
(LO) Signatures of these so called FFLO phases have 
been seen in an atomic gas in one spatial dimension Q, 
but a decisive experiment that observes the FFLO corre- 
lations has not been carried out yet. In three spatial di- 
mensions a phase with nonzero momentum Cooper pairs 
is predicted to be present in the population imbalanced 
Fermi gas with weak interactions [l(>rL4l|. However, the 
transition temperatures in a weakly-interacting Fermi gas 
are very low and are expected to be out of reach with 
present cooling techniques. 

In the three-dimensional strongly interacting mixture 
of ^Li and '"'K atoms a Lifshitz instability towards a 
phase with nonzero momentum Cooper pairs is present 
at an attainable temperature [l^ [ig, but the precise 
form of this superfluid phase has not yet been deter- 
mined. In this Letter we show that at temperatures be- 
low the Lifshitz point the Fermi gas will form the exotic 
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FIG. 1: Phase diagram for the unitary mixture of ^Li and 
""'K atoms (a), where the cross denotes the point for which we 
calculated the atomic densities in Fig. 3. The grey rectangle 
on the horizontal axis sets the scale for panel (b) , where along 
the line of phase transitions to an inhomogeneous superfluid 
is shown. Also shown are the Fermi wavevectors of the ''Li(t) 
and *''K(4) atoms and their difference, all at T = 0. 



LO phase, which shares the properties of both a solid 
and a superfluid and is therefore a special kind of su- 
persolid. This is a very exciting prospect, because even 
though this new kind of superfluid was first proposed by 
Larkin and Ovchinnikov already in 1964 Q, it has never 
been observed in a three-dimensional fluid. In the ^Li- 
'^''K mixture important steps have already been made 
experimentally. The Feshbach resonances of the mixture 
were known for some time [U [3 > now also expan- 
sion under the influence of strong interactions from an 
optical dipole trap has been realized [19], a large atom 
number dual-species magneto-optical trap was built [20| , 
and the repulsive polaron has been studied in this mix- 
ture [2l|. At present, therefore, an ultra-cold ^Li-^'^K 
mixture at unitarity is to the best of our knowledge the 
only system available in the laboratory to observe such a 
three-dimensional supersolid and study its properties in 
detail with the accuracy of atomic physics. 

The LO state is also expected to be present in a num- 
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ber of other condensed-matter systems, such as neutron 
stars [121 and heavy-fermion systems [2^. In the con- 
text of colorsuperconductivity a large body of theoret- 
ical work has been concerned with different inhomoge- 
neous phases (24h2^] . all considering weak coupling and 
using Ginzburg-Landau expansions of the free energy. 
Although continuous (second-order) instabilities can be 
determined exactly using this expansion, the downfall is 
that ruling out discontinuous (first-order) transitions is 
impossible and one is moreover forced to remain close to 
the normal-supcrfluid transition. 

In this Letter we calculate for different inhomogeneous 
superfluids the full inhomogeneous thermodynamic po- 
tential within mean-field theory, with which both con- 
tinuous and discontinuous phase transitions can be de- 
scribed. The Lifshitz point found using this mean-field 
theory remains present when adding screening and selfen- 
ergy effects and its position only changes quantitatively 
[Tsl [T6| . Adding these effects for the mass-balanced case 
reproduces the experimentally found phase diagram [l^ , 
from which it can be concluded that for these purposes 
the mean-field calculation already contains the relevant 
physics, even at unitarity. 

We determine the phase diagram of the mixture of ^Li- 
'^^K atoms for a majority of atoms, see Fig[lja), and 
especially what the structure of the inhomogeneous su- 
perfluid is. For a majority of ^Li atoms it was shown al- 
ready that a tricritical point is present, below which a dis- 
continuous phase transition towards a homogeneous su- 
perfluid, and thus phase separation, occurs [l^- We find 
for small negative polarizations P = (n-j- — n^)/ (% — jt-j,) , 
where f^-|-(4,) denotes the density of ^LiC'^K) atoms, a con- 
tinuous transition from the normal state (N) to a homo- 
geneous superfluid (HSF), up to the Lifshitz point (LP). 
For larger majorities of atoms there is a continuous 
phase transition to an inhomogeneous superfluid phase 
where translational symmetry is broken in one spatial di- 
rection (LO^) and for even larger population imbalances 
the continuous transition is to a phase where the transla- 
tional symmetry is broken in three directions (LO^). In 
FiglIJa) the temperature T is scaled by the critical tem- 
perature Tcfl a.t P — and the polarization P is scaled 
by the critical polarization \Pc,o\ at T = 0, because in- 
teraction effects in first instance only shift the location 
of the Lifshitz point [l^, [IB] • 

In addition, after explaining how we obtained the 
above phase diagram, we determine the densities for the 
two atomic species in the LO^ phase at the conditions of 
the cross in FiglTfa), which are now position dependent 
with a modulation that can be as big as 10% over one pe- 
riod. This is therefore a convenient signature for the LO^ 
phase in experiment, but also other proposals to observe 
this exotic state of matter have been made [H, Is^l • 

Thermodynamic Potential. — We calculate the grand- 
canonical thermodynamic potential f2 for the two- 
component Fermi gas with mass and population imbal- 
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where A(x, r) is the complex pairing field and ^^(x, r) 
are fermionic fields with mo- the mass and /lo- the chemical 
potential of the fermionic atoms. The last two terms in 
the above action represent respectively the annihilation 
and creation of a Cooper pair, here consisting of a ^ Li and 
a atom. The microscopic interaction strength Vq in 
the above action can be replaced in favor of the effective 
interaction proportional to the s-wave scattering length 
a [Tg'I, which can be controlled using a Feshbach reso- 
nance. In this Letter we always consider unitarity, where 
\a\ = oo. From the action the partition function can be 
calculated, Z = / d[0*]d[0]d[A*]d[A] exp(-S'[(/., A]/?i), 
and subsequently the thermodynamic potential can be 
obtained via $7 — —ksT log Z , where ksT is the thermal 
energy with fc^ Boltzmann's constant. Densities can in 
turn be calculated from the thermodynamic potential via 
Ua- = —dQ/dfia-. In the partition function Z the path in- 
tegral over the fermionic fields 4> can be evaluated, after 
which only the path integral over the pairing field A is left 
with in the integrand an exponent of an effective action 
S"^^[A]. In the mean- field approximation the pairing field 
is now replaced by its most probable value and the re- 
maining path integral omitted. In practice the thermody- 
namic potential is then obtained by making an ansatz for 
the pairing field and thus for the superfluid phase. The 
thermodynamic potential describing the phase transition 
from a normal gas to a homogeneous superfluid is calcu- 
lated by taking for the pairing field A(x, r) = Aq. Fulde 
and Ferrel considered a plane wave, A(x, r) — Aqc^'^'^, 
in which case the Cooper pairs have nonzero momentum 
q but the superfluid density is still position independent 
and thus it is a homogeneous phase. Larkin and Ovchin- 
nikov assumed a standing wave, A(x, r) = Aq cos(q • x). 
In this case not only the Cooper-pair phase but also the 
superfluid density becomes position dependent, which 
turns out to be always more favorable than the plane 
wave ansatz. In the latter case the superfluid is truly 
inhomogeneous. Obviously, many other choices for the 
Cooper-pair wavefunction can be made, depending on 
the situation of interest. 

In the case of a position-dependent pairing field the 
thermodynamic potential not only depends on the am- 
plitude Aq but also on the lattice wavevector |q| = q oi 
the order parameter, ri(Ao, q). For all ansatzes we made 
for the order parameter we find at the same chemical 
potential and temperature a continuous phase transition 
to the superfluid phase. Only after comparison of the 
different thermodynamic potentials for different pairing 
fields the most favorable configuration of the superfluid 
is found. In Fig|2] we show the thermodynamic poten- 
tials for different Cooper-pair ansatzes along the line of 
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FIG. 2: The thermodynamic potential along the normal to 
superfluid transition hne from Fig[lja). At first the LO^ 
(solid line) gives the lowest thermodynamic potential, while 
for larger polarizations this is true for the LO^ phase (dashed 
line) . An LO'^ phase (dotted) or more complicated configura- 
tions (dashed-dotted and dashed-double-dotted) never form 
the equilibrium state of this system. 



phase transitions to the inhomogeneous superfluid phase. 
It can be seen that, inside the inhomogeneous superfluid, 
first the phase with only one broken translational sym- 
metry, the LO^ ansatz, gives the lovifest thermodynamic 
potential. Further along the line of phase transitions it 
becomes more favorable to break translational symmetry 
in three spatial dimensions (LO^), corresponding to the 
following ansatz for the order parameter 

A3(x) = Ao [cos{qx) + cos(gj/) + cos(qz)] . (2) 

Other phases are never favorable. By comparing ther- 
modynamic potentials for different lattice structures we 
completed the phase diagram for the Fermi mixture of 
^Li and atoms, see Fig Ufa). The lattice wavevector, 
i.e., the value of q for which the thermodynamic poten- 
tial f2(Ao,q) has a global minimum, along the line of 
the phase transition is shown in Figlljb). It can be seen 
that for T = it differs from the difference in Fermi 
wavevectors, which is due to the fact that we are in the 
strongly-interacting regime. We found that inside the 
LO'^ superfluid phase the global minimum of the ther- 
modynamic potential shifts continuously from a nonzero 
to a zero lattice wavevector q, in which case the superfluid 
is homogeneous again. This criterion for the transition 
between the inhomogeneous and homogeneous superfluid 
is equivalent to the vanishing of the domain- wall energy 

mm- 

LO^ Green's Function. — We found that for one point 
on the normal to superfluid transition line the thermo- 
dynamic potentials for the LO^ and the LO'^ phase are 
equal, see FigUfa) and FigH] Within the superfluid re- 
gion the transition line between these two different inho- 
mogeneous phases can be obtained by comparing thermo- 
dynamic potentials, but now that we know which phases 
to compare we can also calculate the transition line more 



directly. In order to do so, we split the order parameter 
for the LO^ phase 



A(x) = Alo(-z) + Ai[cos(ga;) + cos{qy)] 



(3) 



where Ai^o{z) = Aocos(qz). We now expand the ther- 
modynamic potential for the cubic phase in the perpen- 
dicular order parameter A^ 



^ n{Ao,q) + a±{Ao,q)\A±\ 



(4) 



At the moment a± changes sign the minimum of fia shifts 
to a non-zero value of A^, which means the LO'^ phase is 
more favorable than the LO^ phase. In order to calculate 
a± we need to determine the Green's function for the 
LO^ phase, since 



V_ 



ksT 
2|A,P 



Tr[GLoAiGLoAj 



(5) 



where V is the volume and A±{x,y) — (Aj^/?i)[cos(qx) + 
cos{qy)]ax, with ax the first Pauli matrix in Nambu 
space. The trace in the above equation is taken over 
real space, imaginary time and Nambu space. The in- 
verse Green's function Glo^^ is known from Eq.(IT]) and 
can be expanded in its energy eigenmodes, which are the 
Bogoliubov quasiparticle wave functions. In this man- 
ner we obtain the Green's function in the LO^ phase, 
which contains the inhomogeneities of this phase in an 
exact way. The quasiparticle wavefunctions for the 
Larkin-Ovchinnikov phase we calculate by solving the 
Bogoliubov-de Gennes equation 




huJri 



M„(x) 

w«(x) 
(6) 



where u and v are the quasiparticle coherence factors 
and n is the band index labeling the energy eigenmode. 
The inverse Green's function, on the left-hand side of the 
equation, contains the LO^ order parameter. Since this 
pairing field has a periodicity the quasiparticle wavefunc- 
tions are according to Bloch's theorem of the form 



u„(x) 

Vn (x) 



Wn,k(z) 



(7) 



where Un^k{z) and Vn,k{z) are periodic functions with the 
same periodicity as the gap Ai,o{z) and are thus only 
periodic in the z-direction, which is reflected also in the 
single band index n. 

We use the equilibrium state of the thermodynamic 
potential for the LO^ phase as input to calculate a± in 
order to determine the transition line between the LO^ 
and LO'^ phase, see Fig[T]Ja). 

Inhomogeneous Densities. — From solving the 
Bogoliubov-de Gennes equation in Eq.® we know 
the quasiparticle coherence factors u and v for the 
LO^ phase. With these wavefunctions it is possible to 
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FIG. 3: Densities for tiie ''Li atoms (a) and the ''"K atoms 
(b) scaled by the total density n — ^^(0) + ri^(O) in the LO^ 
phase as a function of position for the point denoted with a 
cross in Fig[IJa). 

calculate the particle densities for the ^Li and the 
atoms as a function of position. Namely, the densities 
are given by 

n.k 

+<kW[l-A^-^(?i^n,k)]}. (8) 

We take the vector k to lie in the first Brioullin zone in 
momentum space, which is here infinitely large in the x 
and y direction and has size 2q in the z direction. In FiglS] 
the inhomogeneous densities in the LO^ phase are shown 
for the point denoted with a cross in the phase diagram 
FigUJa) and it can be seen that the density of the minor- 
ity atoms follows the form of the Cooper-pair ansatz. In 
contrast, the density of the atoms has its maximum 
at the position where the Cooper pair has a minimum. In 
this way the system relaxes the frustration caused by the 
population imbalance, instead of phase separating, which 
occurs for a majority of ^Li atoms [l^, We find that 
the modulation in the density of the ^Li atoms is about 
10%, as can be seen in FiglSja), whereas the modula- 
tion for the majority atoms is much smaller. The density 
modulations can be visible in an experiment, for exam- 
ple using Bragg spectroscopy. A homogeneous superfluid 
does not show Bragg peaks and thus the appearance of 
any Bragg peak is a sign of inhomogeneities in the super- 
fluid phase. 

Discussion and Conlusion. — For the mixture of reso- 
nantly interacting ^Li and ^'^K atoms we completed the 
phase diagram, which contains inhomogeneous superfluid 
phases. These superfluid phases have a crystalline order, 
just as a solid, and are thus special kinds of supersolids. 



Usually a supersolid is defined as a phase with both long- 
range diagonal and off-diagonal order in the one-particle 
density matrix. The LO^ and LO'^ only have diagonal 
order in the one-particle density matrix but there is both 
long-range diagonal and off-diagonal order in the two- 
particle density matrix. This means that the inhomo- 
geneous superfluid is a supersolid phase, but not in the 
usual sense. 

The results presented in this Letter were calculated 
making use of two approximations. The first one is the 
mean-field approximation, which is known to give good 
qualitative results. Although our results are only reason- 
able estimates quantitatively, we have great confidence 
that, as in the mass-balanced case, the phase diagram 
does not change qualitatively if screening and selfenergy 
effects are included. However, fluctuations can have a 
great effect on the density modulations. In order to 
observe a spatial periodicity in the densities true long- 
range order is most desirable. But due to fiuctuation 
effects only alg ebraic long-range order will exist in the 
LO^ phase |l4l [3ll|. making it more difficult to observe 
the periodicity. 

The second approximation we make is that we take into 
account only a single magnitude of the lattice wavevector 
q. For the continuous transition from the normal gas to 
the superfluid state this is exact, whereas it is an approx- 
imation inside the superfluid phase. Deep in the super- 
fluid phase the pairing held is expected to take a more 
complicated form and to depend on a range of wave vec- 
tors • If in the phase diagram the chemical-potential 
difference, instead of the polarization, is on the x-axis all 
transitions take place close to the normal to superfluid 
transition. Therefore, it is not expected that taking into 
account a more complicated form for the gap will change 
the phase diagram qualitatively. Moreover, our calcula- 
tion is variational and thus gives a conservative estimate 
of the supersolid region of the phase diagram. 

In this Letter we presented results for a homogeneous 
system and we did not consider a conflning trap. To take 
into account the trap for the atoms, a local-density ap- 
proximation can be made. The trap felt by the atoms 
depends on the species of the atoms, which means that 
for two atom species two different traps have to be taken 
into account '19']. The trap can have an advantageous 
influence. Namely, if the periodicity of the Cooper pair 
wavefunction is pinned by the shape of the trap it be- 
comes easier to observe density modulations. 

This work is supported by the Stichting voor Funda- 
menteel Ondcrzoek der Materie (FOM), the Nederlandse 
Organisatie voor Wetenschaplijk Onderzoek (NWO). 
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